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We present a method to non-perturbatively determine the parameters of the on- 
sheh, 0(a)-improved relativistic heavy quark action. These three parameters, mo, Ci 
and cb = ce are obtained by matching finite-volume, heavy-heavy and heavy-light 
meson masses to the exact relativistic spectrum through a finite- volume, step-scaling 
recursion procedure. We demonstrate that accuracy on the level of a few percent can 
be achieved by carrying out this matching on a pair of lattices with equal physical 
spatial volumes but quite different lattice spacings. A fine lattice with inverse lattice 
spacing 1/a = 5.4 GeV and 24^ x 48 sites and a coarse, 1/a = 3.6 GeV, 16^ x 32 
lattice are used together with a heavy quark mass m approximately that of the 
charm quark. This approach is unable to determine the initially expected, four 
heavy-quark parameters: mo, Q, cb and ce- This apparent non-uniqueness of these 
four parameters motivated the analytic result, presented in a companion paper, that 
this set is redundant and that the restriction ce = cb is permitted through order 
a\p\ and to all orders in ma where p is the heavy quark's spatial momenta. 

PACS numbers: 11.15.Ha,12.38.Gc,12.38.Lg,14.40.-n 



I. INTRODUCTION 

The study of flavor physics and CP violation plays a central role in particle physics. In 
particular, many of the parameters of the Standard Model can be constrained by measure- 
ments of the properties of hadrons containing heavy quarks. However, to do this one needs 
theoretical determination of strong-interaction masses and matrix elements to connect the 
experimental measurements with those fundamental quantities. Lattice quantum chromo- 
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dynamics (QCD) provides a first-principles method for the computation of these hadronic 
masses and matrix elements. However, lattice calculations with heavy quarks present special 
difficulties since in full QCD calculations, which properly include the effects of dynamical 
quarks, it is impractical to work with lattice spacings sufficiently small that errors on the 
order of ma can be controlled. These problems are addressed by using a number of improved 
heavy quark actions designed to control or avoid these potentially important finite lattice 
spacing errors. The results of recent calculations of basic parameters of the Standard Model 
can be found in the lattice heavy quark reviews of Refs. 

A variety of fermion actions have been used in lattice calculations involving heavy 
quarks . These include heavy quark effective theory (HQET) (for which 

the static approximation is the leading term) and non-relativistic QCD (NRQCD). These 
methods have significant limitations: NRQCD has no continuum limit, and although HQET 
has a continuum limit, it cannot be applied to quarkonia. While systems involving bottom 
quarks may permit a successful expansion in inverse powers of m, this is likely not true for 
systems including a charm quark. 

A third approach, the one adopted here, is the Fermilab or relativistic heavy quark (RHQ) 
method in which extra axis-interchange asymmetric terms are added to the usual 

relativistic action. As is discussed below, this action can accurately describe heavy quark 
systems provided the improvement coefficients it contains are properly adjusted. As the 
heavy quark mass decreases, this action goes over smoothly to the order a improved fermion 
action of Sheikholeslami and Wohlert (SW) Thus, it seems appropriate to refer to this 
as the relativistic heavy quark method since it retains the relativistic form of the Wilson 
fermion action (with the exception that lattice axis interchange symmetry is broken) and 
approaches the standard relativistic action as ma becomes small. 

nn 

As is discussed in the original papers [12, Il3| and considered in detail in our companion 
paper this approach builds upon the original work of Sheikholeslami and Wohlert, 
extending it to to the case of a possibly very heavy quark with mass m > 1/a but restricted 
to a reference system in which these quarks are nearly at rest. Such a situation can be 
described by a Symanzik effective Lagrangian which contains terms of higher dimension 
than four which explicitly reproduce the finite lattice spacing errors implied by the lattice 
Lagrangian. 

In this RHQ approach, the continuum effective Lagrangian is imagined to reproduce 
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errors of first order in a\p\ and all orders in ma or p^a where {p,p^) is the heavy quark 
four-momentum. Such an effective action will contain many terms, including those with 
arbitrarily large powers of the combination aD^ where the gauge-covariant time derivative 
will introduce a factor of p° and so cannot be neglected. As described in Ref. [12] and 
discussed in detail in our companion paper, this Lagrangian can be greatly simplified by 
performing field transformations within the path integral for the effective theory. These 
field transformation do not change the particle masses predicted by the theory and, as is 
discussed in Ref. they effect on-shell spinor Green's functions only by a simple, Lorentz 
non-covariant, 4x4 spinor rotation. (The use of the equations of motion in Ref. is 
equivalent to the above field transformation approach.) 

As is shown in our companion paper after this simplification, the resulting effective 
Lagrangian contains only three parameters: the quark mass m'^, an asymmetry parameter 
describing the ratio of the coefficients of the spatial and temporal derivative and a general- 
ization of the Sheikholeslami and Wohlert, csw fo the case of non-zero mass which we refer 
to as Cp. Here the superscript c indicates that these are the coefficients that appear in the 
continuum effective Lagrangian. If these three, mass- dependent parameters can be tuned to 
physical continuum values (m^, 1 and respectively) by the proper choice of mass- dependent 
coefficients in the lattice action, then the hadronic masses computed in the resulting theory 
will contain errors no larger than (pa)^. 

In addition, a new parameter S multiplying a non-covariant 7 ■ p in the 4x4 spinor matrix 
mentioned above will be needed to realize truly covariant on-shell Greens functions. Here 5 
will depend on the (usually composite) fermion operator being used even for a fixed action. 



As is discussed below and in detail in Ref. [15|], this is one fewer parameter than found in 
the previous work of the Fermilab and Tsukuba groups. 

Thus, in our calculation we use the relativistic heavy quark lattice action: 

-^lat = ^n' (70^0 + C7 • ^ + "^0 - 2^t^0 - 2^sD 

n' 



4 ^ ■ ^ 2 



where 



D^i^{x) = i [f/^(x)^(x + fi)- Ul{x - fi)ij{x - fi)] (2) 
Dl^Pix) = \u^ix)^lj{x + + Ulix - fi)^{x -fi)- 2^(x)l (3) 
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ss' [Usf,{x)Us'u{x + sfi) U^sf,{x + sfi + s'O) 
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s,s'=±l 



(4) 



and the 7 matrices satisfy: 7^ = 7^, {7/^,7i/} = 25^i, and cr^j/ = |[7'^,7'^]. Written in 
this standard form, there are six possible parameters that can be adjusted to improve the 
resulting long-distance theory, three more than are needed. We begin by making the choice 
Ts = C and Tt = 1. This leaves four parameters whose effects we can study. In the following 
we will investigate the non-perturbative effects of these four parameters, mo, cb and ce- 
However, when determining an improved RHQ action non-perturbatively, we will impose 
the further restriction cb = ce, making the improved action uniquely defined at our order 
of approximation. 

The different coefficient choices of improved lattice action by the Fermilab and Tsukuba 
groups yield two distinct sets of coefficients for the action. These are summarized in Table HI 
The coefficients in each approach have been calculated by applying lattice perturbation the- 
ory at the 0{a)-improved, one-loop level to the quark propagator and quark-quark scattering 
amplitude jl^ III | • 

In this paper, we will propose and demonstrate a non-perturbative method for deter- 
mining these coefficients based on a step-scaling approach, which eliminates all errors of 
0((7^"). Step scaling has been used in the past to connect the lattice spacing accessible 
in large-volume simulations with a lattice scale sufficiently small that perturbation theory 
becomes accurate Non-perturbative matching conditions are imposed to connect the 
original calculation with one performed at a smaller lattice spacing a[ = ae. Iterative match- 
ing of this sort with n steps then connects the theory of interest and a target lattice theory 
defined with lattice spacing a'^ = ae". This may require a number of steps n which is not too 
large since while the coupling constant decreases only logarithmically with the energy scale, 
that energy scale increases exponentially with n. For example, if a final comparison with 
order g"^ perturbation theory is used, we expect an error of order {g')'^ ~ In(a^)^^ ~ 
where g' is the bare coupling for the finest lattice. 

The situation for heavy quarks is even more favorable. Here the target, comparison theory 
does not need to have such small lattice spacing that perturbation theory is accurate. In 
fact, this theory can be treated non-perturbatively provided the final lattice spacing a'^ is 
sufficiently small that simulations with an ordinary 0(a)-improved relativistic action will 



give accurate results This implies that the size of the error will be of order {ma'Y ~ 
or n ~ (error) In the work reported here we will match the step-scaled heavy quark 
theory with an 0(a)-accurate lattice calculation performed using domain- wall fermions. 

A critical question in developing such a step-scaling approach is to decide upon the actual 
quantities that will be "matched" when comparing two theories defined with different lattice 
spacings but which are intended to be physically equivalent. Among the quantities which 
might be matched are the Schrodinger functional 'iS', off-shell Green's functions defined in 



the RI/MOM scheme 20| or physical masses and matrix elements at finite volume. 



Our first approach to this topic was to investigate the off-shell RI/MOM scheme since this 
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and also 



method had worked well in earlier light-quark calculations, see e.g. Refs. 
permits a direct comparison with quantities defined in perturbation theory. We were able 
to define RI/MOM kinematics which lay within the regime of validity of the effective heavy 
quark theor y d escribed above and to carry out a tree-level calculation of the amplitudes of 
interest j23l . l24| . However, the increased number of parameters needed in the effective theory 
to describe off-shell amplitudes, the need to work with gauge-non-invariant quantities and 
the difficulty of computing "disconnected" gluon correlation functions ultimately made this 
approach appear impractical. 

In this paper, we adopt the third method mentioned above and determine the coefficients 
in the heavy quark effective action appearing in our step-scaling procedure by requiring that 
the physical, momentum-dependent mass spectrum of two physically equivalent theories 
agree when compared on the same physical volume. Since the step-scaling approach requires 
physically small volumes be studied, these spectra will be significantly distorted by the effects 
of finite volume and it is important that these effects be the same in each of the theories being 
compared — thus the need to compare on identical physical volumes. By comparing more 
physical, finite-volume quantities (as many as seven) than there are parameters to adjust 
(three), we also have an over-all consistency test of the method. Finally, as described above, 
at the smallest lattice spacing, we compute the quantities being compared using a standard 
domain wall fermion action which is has no order ma errors and accurate chiral symmetry. 
We assume that at this smallest lattice spacing the explicit errors of order (ma)^, present 
in the domain wall fermion calculation, are sufficiently small to be neglected. Preliminary 
results using this method were published in Ref. j^]. 

The structure of this paper is as follows. We introduce our on-shell approach to determine 
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the coefficients in the relativistic heavy quark action via step-scahng both in the quenched 
approximation and for full QCD in Sec. |nj In this paper we will work in the quenched 
approximation and explicitly carry out the first stage of matching between a fine and a 
coarse lattice in order to determine the feasibility of this approach and the accuracy that 
can be achieved. Specifically the "fine" lattice uses 1/a = 5.4 GeV and the "coarse" lattice 
1/a = 3.6 GeV. (A second matching step, evaluating the first coarse lattice action on a 
larger physical volume and matching to an even larger lattice spacing, 1/a = 2.4 GeV, 
is now underway.) Section ITTTl lists the parameters used in this calculation, describes our 
determination of the lattice spacing and discusses our method for obtaining the physical 
heavy-heavy and heavy-light spectrum. 

The problem of determining the parameters to be used in the coarse-lattice effective 
action which will reproduce the fine lattice mass spectrum is studied in Sec. IIVI and the 
dependence of this spectrum on the four parameters mo, and ce presented. We are 

unable to determine these four parameters with any reasonably precision, a conclusion we 
now understand since only three parameters are required to determine the mass spectrum 
to order a\p\ and all orders in (ma)" [15]. We then restrict the parameter space to = 
Ce, as is justified theoretically, and find that the resulting three parameters can now be 
determined quite accurately. In Sec. |V|we compare our result with both perturbative and 
non-perturbative determinations of the quark mass and the one-loop lattice perturbation 
calculation of the lattice parameters (, cb and ce performed by Nobes 2^. Section. IVII 
presents a summary and outlook for this approach. 



II. STRATEGY 



We propose to determine the three coefficients mo, C and cp = cb = ce in the RHQ lattice 
action of Eq. Qby carrying out a series of matching steps. We begin with a sufficiently fine 
lattice spacing that no heavy quark improvements are needed (ma <^ 1) and a conventional 
light-quark calculation will give accurate results. We then carry out a series of calculations 
using the RHQ lattice action of Eq. ^on lattices with increasingly large lattice spacing and 
increasingly large physical volume. When we increase the lattice spacing at fixed physical 
volume, we perform calculations at both lattice spacings on identical physical volumes and 
require that the resulting finite- volume, heavy-heavy and heavy-light energy spectrum agree 
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when these particles are at rest or have small spatial momenta. When we increase the lattice 
volume at fixed lattice spacing, we simply use the parameters, previously determined, in a 
calculation now on the larger volume. The first same-physical-volume matching of the energy 
spectrum is done between the heavy quark theory and a conventional fine-lattice calculation 
done with domain wall fermions. An example of this finite volume, step-scaling recursion is 
shown in Fig. ^ 

Specifically, we will calculate the pseudo-scalar (PS), vector (V), scalar (S), and axial- 
vector (AV) meson masses in the heavy-heavy (hh) system and pseudo-scalar and vector 
masses for the heavy-light (hi) system. We will work with the following mass combinations: 

• Spin-average: m'^!^ = \ {m^ + Smy'*^ , mJ^l = j (rripg + 3my ^ 

• Hyperfine splitting: m(jg = rriy'^ — mp|, m^^ = my — mpg 

• Spin-orbit average and splitting: m^^^ = \ (m^'^ + 3m^y^ , m^l^^ = m^y — mg'' 

• Mass ratio: mi/m2 where E"^ = ml + ^p^, with mi the rest mass and m2 the kinetic 
mass. 

By examining these seven quantities we should be able to determine the three parameters 
mo, C ^iid cp and also check the size of the scaling violations. 

The first step in this program calculates these seven quantities using the domain wall 
fermion action on the fine, 24^ x 64 lattice with 1/a = 5.4 GeV (I). Next, these seven 
quantities are computed a second time using a coarse, 16'^ x 32 lattice with 1/a = 3.6 GeV 
and, therefore, the same physical volume. This is calculation II. The three heavy-quark 
parameters entering this coarse-lattice calculation must be adjusted so that these seven 
quantities agree between calculations I and II. It is these calculations that are carried out 
in this paper using the parameters given in Table |nl 

Third, we expand the volume of calculation II to 24'^ x 48, while keeping all other parame- 
ters fixed. The results of this third, expanded volume calculation (III) can then be matched 
with a fourth calculation which has a lattice spacing larger by a factor of 3/2 (IV). The 
simulation parameters for this second matching step are given in Table IIIII By repeating 
this pattern, we can extend the calculation to quenched lattices with the desired volume 
where serious, infinite- volume, charm physics may be studied. 
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In this paper, we demonstrate only the matching between calculations I and II. The lead- 
ing heavy quark discretization error in calculation I is (ma)^ ^ 4%, making it the dominant 
systematic error on the fine lattice result. Of course, this error can be reduced in future 
calculations by choosing a fine lattice that has an even smaller lattice spacing and corre- 
spondingly smaller physical volume. Without improvement beyond the usual Sheikholeslami 
and Wohlert term, the leading heavy quark discretization error on the coarse lattice is ex- 
pected to be (ma)^ ^ 10%. However, once we introduce the improved lattice action of Eq. ^ 
and properly tune the coefficients, we should be able to reduce the error to (a|p|)^ ~ 1%. 

As will be demonstrated in the remainder of this paper, this proposed step-scaling method 
works well and offers a feasible approach to heavy quark calculations with accurately con- 
trolled finite lattice spacing errors. However, unless we can move beyond the quenched 
approximation used here, this method will be of only limited utility. Using this method for 
full QCD will, of course, be more computationally demanding because each of the two sets 
of lattice configurations generated for this matching process must be obtained from a full, 
dynamical simulation including the quark determinant. However, such an approach could 
become prohibitively expensive if the value of the light dynamical quark masses, mughta, 
must to decrease toward zero with decreasing a. 

Fortunately, such small dynamical quark masses are not required in this step-scaling 
approach. The combination of the usual gauge and light-quark action together with the 
effective lattice action of Eq. Q defines a complete physical theory, including the properties 
of heavy quarks, that is unambiguously specified at short distances A with a <^ A ^ l/mught. 
Recall that in the continuum such a local field theory is typically defined in a mass- and 
volume-independent fashion. Short-distance renormalization conditions are imposed to fix 
the theory in a manner that is insensitive to quark masses and space-time volumes. Similarly 
our fine-lattice theory, viewed as a function of the bare input lattice parameters also defines 
such a mass- and volume-independent theory. 

Given sufficient computer power, the implications of this theory could be worked out 
on arbitrarily large spatial volumes and for arbitrarily small masses. The results would be 
well-defined functions of the bare input parameters which would require no adjustment as 
the quark mass and spatial volume were varied. The lattice spacing could be determined in 
physical units by comparing to Aqcd as determined from a vertex function at short distances 
with the light quark masses having a negligible effect. 



9 



Replacing the standard light quark action appropriate for our finest lattice with the 
improved lattice action of Eq. ^does not change this situation. The parameters in the heavy 
quark action could be evaluated or renormalized by examining Green's functions evaluated 
at non-exceptional momenta without infra-red or light-quark mass sensitivity [s^l^]. We 
would still be working with a short-distance-defined field theory that will give meaningful 
predictions as a function of lattice volume and light quark mass. Thus, when comparing 
two such effective theories defined at two different lattice spacings we are free to use any 
lattice size L and dynamical quark mass miight we find convenient provided L ^ a and 
'flight ^ 1/a. In fact, if mught is sufficiently small that it does not effect the finite- volume 
heavy-quark spectrum being compared, we need not even use the same quark mass in the 
two calculations being compared! Of course, this is likely a quark mass that is expensively 
small and a better strategy is to work with sufficiently small spatial volume and sufficiently 
heavy dynamical quark mass that they do effect the quantities being matched and must be 
given equivalent physical values in each of the calculations being compared. 

We conclude that employing the procedure developed here in full QCD, while difficult, is 
practical and well within the reach of present computer resources. Just as our step-scaled 
lattice spacing decreases and we move to increasingly smaller spatial volumes, we should also 
move to increasingly heavier quark masses. In both cases finite volume and finite dynamical 
quark mass effects are distorting the spectrum being compared, but these distortions are 
entirely physical and must be accurately described by the effective actions being compared. 



III. SIMULATIONS 



We performed this calculation on a 512-node partition of the QCDOC machine located 
at Columbia University. We used the Wilson gauge action since for this case the relation 
between lattice coupling and lattice spacing has been thoroughly studied 

The gauge configurations were generated using the heatbath method of Creutz js^, 
adapted for SU{3) using the two-subgroup technique of Cabibbo and Marinari jsil. The 
first 20,000 sweeps were discarded for thermalization and configurations thereafter were 
saved and analyzed every 10,000 sweeps. We examined the auto-correlation between config- 
urations for both the standard 4-link plaquette and the hadron propagators evaluated at a 
time separation of 12 lattice units. Here we use the standard autocorrelation function p{t) 
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as 



1 



Ntot-t 



Pit) 



Ntot - 1 



E (0(j)-0)(0(t + j)-0) 



(5) 



i 



and identify the autocorrelation time as the size of the region near t = in which p{t) ^ 0. 

For the case of the plaquette (which was calculated every sweep) we found an auto- 
correlation time of approximately 3 sweeps. We studied the propagator correlations using 
two 24^ X 32, (3 = 6.638 test calculations: in the first, the propagators were computed on 
120 configurations, separated by 5 sweeps and in the second on 40 configurations separated 
by 50 sweeps. The resulting correlations functions for five different hadron propagators 
evaluated at a temporal source-sink separation of twelve lattice units are shown in Fig. El 
Auto-correlation on the scale of 100 sweeps can be seen in the data sampled every five 
sweeps. Essentially no autocorrelation can be seen for the propagators sampled every fifty 
sweeps. Since we used a total of 100 such lattice configurations, separated by 10,000 sweeps, 
for all of the quantities discussed in this paper, we will assume that quantities calculated on 
different configurations are statistically independent. 



Four different quantities with a meaningful physical scale enter each of the two lattice 
calculations that must be matched at a given stage in our step-scaling procedure. Most 
familiar is the distance scale determined by the static quark potential or the chiral limit 
of the light hadron spectrum. This is an important physical scale that will influence even 
small-volume, heavy quark results. As is conventional, we will refer to this quantity as the 
"lattice spacing" . We will find it convenient to determine this from direct calculation of the 
static quark potential. The other three scales are the lattice volume, and the masses of the 
light and heavy quarks. (In our discussion of heavy-light systems we will ignore the strange 
quark and work with degenerate up and down quarks.) 

Of course, the lattice spacing, expressed in physical units, is also important since it gives 
us a direct idea of the size of the discretization errors which we are trying to control. For this 
purpose we don't need great precision (something that cannot be achieved in a quenched 
calculation under any circumstances). We will determine the lattice spacing in physical units 
from the static quark potential evaluated at an intermediate distance to yield the Sommer 
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Lattice scale 
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scale j28| which is then determined from a phenomenological, static quark potential model. 
While the later is not precisely defined, this method has the advantage that it uses only 
pure gauge theory without any fermion action being involved. We could get a physical value 
for the lattice scale from the pion decay constant /,r or the rho meson mass but these are 
more difficult to calculate and may be more sensitive to finite volume and other systematic 
errors. 

Our strategy for choosing the parameters to be used on the fine lattice and then finding 
physically equivalent parameters to be used on the coarser lattice proceeds as follows. We 
first decide on a target ratio of lattice scales determined by the ratio of lattice volumes that 
we intend to use. In the present case a ratio of 3/2 is implied by our choice of 24^ and 16^ 
lattice volumes. Second, we choose the bare lattice coupling on the fine lattice to insure that 
the fine lattice spacing is sufficiently small (here chosen to be 1/a = 5.4 GeV). The stronger 
coupling must then give a lattice spacing larger by a factor of 3/2 than the fine value or 
1/a = 3.6 GeV. While for a dynamical QCD calculation, this would require considerable 
numerical exploration, for a quenched calculation with the Wilson gauge action, we can 
simply refer to extensive earlier work. 

Next we choose the light quark mass to be used on the fine lattice as sufficiently light 
that the heavy-light mesons being studied will be involve a different momentum scale than 
do the heavy-heavy mesons but not so light as to unreasonably increase the computational 
cost. For the calculations reported here we used the domain wall formalism for the light 
quarks and chose the mass nifa = 0.02, one-tenth of the 0.2 heavy quark mass. The light 
quark mass to be used on the coarse lattice is determined by requiring that the light-light 
pseudo-scalar meson have a mass 3/2 times larger than that found on the fine lattice when 
measured in lattice units. 

Finally the heavy quark mass on the fine lattice is estimated to correspond to the bare 
charm quark mass. While in the present calculation we have used the single value mja = 0.2, 
a complete calculation will likely require one or two more masses so that a final interpo- 
lation/extrapolation can be done to make the physical charmed hadron mass agree with 
experiment. The heavy quark mass on the coarse lattice is one of the three heavy-quark 
parameters whose determination is discussed in the next section. 

Let us now discuss the choice of lattice scale in more detail. The static potential is 
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expected to have the following form 

V{R) 



^-1 



a 



R 



(6) 



where R is the separation between the static quarks. The scale implied by the heavy quark 
potential is often specified using the Sommer parameter tq which is defined by the condition 



R 



.2 dV{R) 
dR 



R=ro 



1.65. 



(7) 



This is appropriate on standard size lattices for bare couplings in the range (3 = < 6.57. 
For weaker couplings, P > 6.57 one uses a second, smaller distance scale, Vc defined by 



R''F{R)\R=r. = 0.65 



(8) 



where ^ = 0.5133(24) i29|. While it may be problematic in a quenched calculation, we 
can attempt to determine ro from a phenomenological potential model, which gives r^^ = 
0.395 GeV. 

Reference j^] gives predictions for the resulting lattice spacing when the coupling /5 of 
Wilson gauge action is in the range 5.7 < P < 6.92: 



ln(a/ro) = -1.6804 - 1.7339(/3 - 6) + 0.7849(/3 - 6)^ - 0.4428(/5 - 6)' 



(9) 



With the help of Eq. IHl we can locate the /3 values needed to achieve the desired cutoff 
scales and fine tune it later as necessary. As our final choices, we have (3 = 6.638 for the 
a^^ = 5.4 GeV lattice and (3 = 6.351 for the a^^ = 3.6 GeV one. 

Since the comparison of lattice scales between our two simulations is fundamental to this 
matching program, we have carried out additional calculations to make sure that the lattice 
spacing is correctly selected. This requires a direct calculation of the static quark potential 
on our lattice configurations. 

Recall that the static quark potential can be extracted from the ratio of Wilson loops: 

{W(r,t)) 



V(r) = log 



(10) 



_{W{r,t + l))_ 

where (...) denotes an average over gauge configurations. In order to improve the signal 
and to extract the potential V{r) from smaller time separations, we smear the gauge links 
in the spatial directions according to Ref. |32l | 



SU(3) 



Uk{n) + 

"-^smear 



(11) 
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where k and / each indicate a spatial direction, Psu(3) is an operator that projects a hnk back 
to an SU(3) special unitary matrix, Csmear is the smearing coefficient (set to 0.5 in our case) 
and the smearing procedure is performed rigmear times. More details regarding the algorithm 
can be found in Ref. 3]. In our calculation we found good results for nsmear = 180 for 
the fine lattice and rismear = 60 for the coarse one. 

While we did determine the two scale standards tq and individually for both of our 
lattice spacings, our lattice volumes are somewhat small to permit a comparison with infinite 
volume results. We therefore also determined the ratio of lattice spacings without using the 
Sommer scale by directly comparing the potentials computed on our two sets of lattice 
configurations using the relation 

Vi{n) = V2{n/X)/X + C' 

where A = 02/01 is the ratio of the two lattice spacings. We first fit the static potential on 
the fine lattice to the form given in Eq. IHl determining the parameters C, a and a. Next we 
scaled the resulting fitted function according to Eq. [T21and adjusted the parameters A and C 
in that equation to obtain the best fit to the static potential measured on the coarse lattice. 
Figure El shows a comparison of the potential determined from (3 = 6.351 configurations and 
a scaled and shifted version of the (3 = 6.638 potential. The agreement is excellent and this 
procedure gives an independent value for the lattice spacing ratio of 1.51(2), which agrees 
with what we wanted. 



B. Domain- Wall Fermions 



We will now briefly describe the domain wall fermion calculations that were used for the 
heavy quark on the finest lattice and the light quarks on both lattices. The domain wall 
Dirac operator can be written as 



(12) 



1 

^ E [(1 - l^.)Ux,^,^x+i,,x' + (1 + l^^Ul,^^6^-f,^^] + (M5 - 4)4,x' (13) 



2 
2 



1 - 75)5s+l,s' + (1 + l->)5s-l,s' - 25s, 
(1 - 75)5s,L,-l'5o,s' + (1 + 75)5s,o5l.-1 



(14) 
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where the fifth-dimension indices s and s' lie in the range < s,s' < Lg — 1, M5 is the 
five-dimensional mass and mj directly couples the two walls at s = and s = Lg — 1. It is 
related to the physical mass of the four- dimensional fermions. 

The M5 parameter is optimized by the choice of M5 1 — merit, where merit is the 
critical value of the mass for the 4-dimensional Wilson fermion action. This quantity has 
been calculated perturbatively up to one-loop level for the Wilson gauge action and either 
the Wilson0 or SW 3 

fermion actions. In the quenched approximation, for Wilson 
fermions and with our choices of gauge coupling, we find m^it = —0.495 at /? = 6.638, and 
'^crit = —0.522 at /? = 6.351. Therefore, we use M5 = 1.5 in the DWF action for both our 
P values. 

The DWF action is 0{a) off-shell improved due to the preservation of chiral symmetry, 
and no further improvement in the action or quark fields is performed. The chiral symmetry 
breaking can be measured from the residual mass, which can be computed from the ratio 



amr 



E.(J5",(x,t)7r(0)) 



(15) 



E.(Jf(x,t)7r(0))' 

provided t ^ a. Here J^g is a pseudoscalar density located at the midpoint of the fifth 
dimension. The residual mass has been thoroughly studied, for example, in Ref.|37|] for 
various values of /?, Lg and M5. Those results suggest that the mres values for each of our 
lattice configurations are much smaller than the 0.00124 value determined at /? = 6.0 with 
lattice volume 16^ x 32 x 16 and M5 = 1.8. This indicates that chiral symmetry breaking 
is small and ignoring the contribution of mres in the heavy quark sector will have an effect 
smaller than 0.5%. 

However, there is a limitation to using large values of m^ with DWF. Recall that there 



are two types o: 
caying states [38 



eigenvectors of the hermitian DWF Dirac operator: propagating and de- 



39|. The former, unphysical states have non-zero 5*'^-dimension momenta 



and large Dirac eigenvalues around 1/a. The "decaying" states are bound to the walls of 
the 5*^^ dimension and are the physical states corresponding to the four-dimensional Dirac 
eigenstates in the continuum limit. The gap between these two types of states is controlled 
by the domain- wall height M5. However, as increases, the eigenvalues of the physical 
states increase while those of the propagating states do not. Thus, we must be careful to 
avoid the situation in which the states with the smallest eigenvalues are dominated by these 
unphysical states. Therefore, a careful check on the lowest eigenvalues for the target m/ 
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being used to simulate the heavy quarks on the fine lattice is needed. Figure HI shows the 
5-dimensional eigenf unction, averaged over 4- dimensional space, l^x.sP as a function of 
the 5**^-dimensional coordinate s for the lowest nineteen eigenvalues with various m^: 0.22, 
0.27, 0.37, 0.47. As can be seen in the figure, these first nineteen eigenfunctions appear to 
be physical states bound to the 4-dimensional wall for the first three mass values. However, 
TTLf = 0.047 is sufficiently large that propagation into the 5-dimension can be clearly seen. 
We conclude that our rrif = 0.2 for the heavy quark is safe, well below the region where 
such unphysical states arise. 



C. Spectrum Measurements 

In order to get good signals for the heavy quark states of interest for relatively small time 
separations, a smeared wavefunction source is used for the heavy quark (but not the light 
quark). Here, we adopt the Coulomb gauge-fixed hydrogen ground-state wavefunction: 

^gndir) = e-^/^° (16) 

as the source of the heavy fermion(s) and use a point source for the light fermion (if any). 
At the sink the two propagators are evaluated at the same point and the resulting gauge 
invariant combination summed over a 3-dimensional plane at fixed time, with a possible 
momentum projection factor. An optimized radius, vq, was chosen in the fashion suggested 
in Ref . [4^ . Table IIVI lists all the local meson operators used in our calculation. 

Figure El shows how the plateau in the effective mass plot improves between a point and 
smeared source. The smeared-source meson plateaus are much better those of the local 
source, even for the scalar and axial-vector mesons. 

To constrain the space-time asymmetry parameter (, we also computed the pseudo-scalar 
meson energy in the heavy-heavy sector for the three lowest on-axis momenta: ^(0,0,0), 
^(0,0, 1), and ^(0, 1, 1), where L is the spatial lattice size. The dispersion relation may 
be expanded in momentum as 

2 

E{p) = m, + ^ + 0{p'). (17) 

As we will see, requiring the ratio of static to kinetic mass, m2/mi = 1 is useful for deter- 
mining the coefficient C,. 
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D. Parameters 

Table lists the fixed parameters used throughout this matching stage for the fine and 
coarse lattices. The heavy quark mass was set to approximate that of the charm mass and 
the light quark mass was chosen ten times smaller. The lattice spacing ratio between these 
two lattices is 1.5. The domain wall fermion parameters used on the fine lattice have been 
carefully studied and we find no unphysical states in the chosen mass range as discussed in 
previous section. Table IVTl shows the hadron mass spectrum computed on the fine lattice. As 
can be seen, mi/m2 = 1.02(2) is consistent with 1, indicating that heavy quark discretization 
effects using domain wall fermions are small. One expects that the light quark mass on the 
coarse lattice should be 0.03 and the data for the light-light spectrum with this choice of 
light quark mass is listed in Table Ivnl As one can see from Table IVTTl the light-light meson 
spectra on the coarse and fine lattices agree when compared in the same units, indicating 
that the light quark mass is well tuned on the coarse lattice. 

A complete list of the parameter sets used for the RHQ action on the coarse lattice 
IS given m Table IVTTTl The first 42 sets of data were initial trials chosen to give good 
coverage in parameter space. In order to perform a more systematic analysis, described in 
Sec. IIVI we also collected a "Cartesian" set (sets #43-#66) chosen close to the desired fine 
lattice measurements. These 24 data sets are centered around set #14. The range of each 
parameter in this Cartesian data set was selected so that within that range the estimated 
difference between a linear and quadratic fit would be less than 5% as expected from an 
examination of the first 42 parameter sets. This yields a region that is close to reproducing 
the target fine data and in which a linear approximation should be good: mo = 0.0328 ±0.1, 
cb = 1.511 ± 0.1, ce = 1.538 ± 0.3 and C = 1-036 ± 0.02, which is shown in Figure El 

The 24-set "Cartesian" data will allow us to calculate the first and second derivatives 
directly from the measurements. Note that we have more measurements than we actually 
need. This provides additional checks on our method and the validity of the scaling of 
physical quantities between the coarse and fine lattices. We expect that the total number 
of data sets that we will use for next step of matching will be dramatically reduced. 

Using the methods described in Sec. IIII CI we have measured the pseudoscalar (PS), 
vector (V), scalar (S) and axial- vector (AV) mesons in the heavy-heavy system, and PS 
and V in the heavy-light system. We use combinations of the masses to try to simplify 
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their dependence on the coefficients of the RHQ action. For the heavy-hght system, we use 
the spin-average and hyperfine sphtting; for the heavy-heavy system, we use these and also 
include the spin-orbit average, spin-orbit splitting and the ratio of mi/m2- The resulting 
values for these quantities for each of the 66 data sets are given in Table IIXI 



IV. ANALYSIS AND RESULTS 



The final step in this matching procedure is to determine the parameters in the heavy 
quark action of Eq. [H {mo, cb,ce,C}: that will yield the seven quantities measured on the 
coarse lattice which agree with those determined on the fine lattice. Of course, this might 
be done by "trial and error" and, as can be seen by scaling the numbers in Table IIX[ data 
set #14 comes very close to such a result. However, to fully understand this step scaling 
method (for example to properly propagate errors) it is important to learn in detail how the 
measured spectra depend on these input parameters. 

As a starting point, we will attempt to use a subset of our parameter space chosen so 
that the resulting coarse lattice hadron masses are well fit by a simple linear dependence on 
the heavy quark parameters: 



F" = A + J-X", 



(18) 



where n labels the parameter set while X and Y are 4-dimensional and 7-dimensional column 
vectors made up of the four input heavy-action parameters and the seven computed masses 
or mass ratios, respectively: 



X 



Cb 
Ce 

{ ^ J 



Y 



m 



m 



m 



hh 
hs 

hi 



m 



hi 
hs 



m 
m 



hh 
soa 

hh 



(19) 



y mi/m2 J 

The quantities A and J are a 7-dimensional column vector and a 7 x 4 matrix which represent 
the constant and linear terms in our linear approximation. (In most of the discussion to 
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follow we will work with all seven measured quantities. However, if this number is decreased, 
the vectors Y, A and the matrix J will shrink appropriately.) 

Given a specific group of N of our data sets, X"'|i<j<Ar, we can determine the quantities 
A and J by minimizing an appropriate fo^' such a linear fit: 

N 
1=1 

Here is a 7 x 7 matrix representing a choice of correlation matrix. In the results that 
follow we will use 

N 

iWcU,' = E((>T^ - ^0 {y^'' - Y^J)) (21) 

i=l 

where (...) represents an average over the 100 jackknife blocks obtained by omitting one of 
the 100 measurements with Y^^ the result for that jackknife block and Y^^ the corresponding 
average. Replacing Wc by the simpler, uncorrected error matrix {W^)d,d' = YliLi ^d,d'(y'd' 
had little effect on the final results where cx^ is the usual squared error on the measured 
quantity Y^. Determining the A and J which minimize Xc is straight-forward because this 
is a quadratic function of these 35 numbers and the minimum can be obtained by solving 
35 linear equations. Typically these 35 equations are quite regular, with a stable solution 
even if only a relative few of our data sets are used. 

The use of linearity to determine the desired matching heavy quark parameters is rea- 
sonable if we are working in a region that is close to the right choice for those parameters. 
Once we have determined the matrix J and vector A, we can solve for the coefficients Xc 
that will yield meson masses equal to those found on the fine lattice, Yjr. Here we add the 
subscripts C and to indicate our estimates for the physical coarse-lattice parameters (X^) 
and the coarse-lattice masses (IV) ? scaled from those calculated using the fine lattice. 

Again we minimize a quantity x^^, similar to that given in Eq.EHl However, the fine-lattice 
correlation matrix, W^, which appears in the equivalent version of Eq. |2niis defined through 
a modified version of Eq. Specifically, the fine-data analogue of Eq. EHis used for all but 
the seventh row and seventh column, which correspond to the quantity mi/m2. Since this 
must be unity in a relativistic calculation (and is one within errors for our DWF results), we 
set {Yj:)q = 1 and the corresponding elements of the correlation (Wc)dj^6,6 = (W^c)6,d^6 = 
for < c? < 6. In order that the resulting correlation matrix be invertible, we arbitrarily set 
(Wc)6,6 = 10^^. This has the effect of constraining the coarse-lattice value of mi/m2 = 1. 
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The resulting minimum is again determine by solving a set of linear equations. That solution 
can be written explicitly as: 

Xc = [j^-W^^-jy^-J^-W:F-{Yr-A). (22) 

Finally, to determine the error on the resulting heavy quark parameters Xq we add in 
quadrature two different sources of error. To compute the first, we use the average value 
of the masses IV, deduced from the fine-lattice calculation, which together with jackknifed 
results for J and A gives us the error on Xq coming from the statistical fluctuations in 
the coarse lattice data. We then estimate the statistical error coming from the fine lattice 
calculation by using the average values for J and A in Eq. |221and the jackknifed values of YJr 
to determine the resulting fluctuations in the resulting heavy quark parameters, Xq caused 
by the statistical errors in the determination of Yjr. 

A. Four-Parameter Action 

As is suggested by the large number of data sets listed in Table I VI III we had greater 
difficulty than expected in determining the four parameters mo, c^, ce and Q. Typically, 
reasonable choices of a subset of the parameter sets from the initial group of 42 parameter 
sets listed in Table IVlnl gave similar values for the final heavy quark parameters. However, 
the derivative matrix was typically quite singular and the resulting parameters, especially 
cei not well determined. In an attempt to make this process more deterministic, we collected 
the 24 Cartesian data sets from which we could determine the matrix of derivatives J from 
simple differences. The result for J agreed very well with that typically determined from 
the fitting procedure describe earlier to the less regular parameter choices in our first 42 
data sets. We conclude that this linear description of our coarse lattice data is a good 
approximation. For simplicity, we present only the results from this final determination of 
J and A from the Cartesian data. 

Specifically, the twenty four parameter choices within our Cartesian data set (#43-7^66) 
use parameters of the form 

X'^ = Xi + a{n)A^■ (23) 

Here the quantity {o"(7^)i}o<i<3 determines the first sixteen parameter choices, where a{n)i = 
^_iynt{n/2^) ^ the expression int(a;) represents the integer part of the number x and the index 
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n varies between and 15. The four parameter increments Aq = 0.1, Ai = 0.1, A2 = 0.3 
and A3 = 0.02 were listed earlier and are displayed in Figure El The remaining eight data 
sets use the values o"(16 + n)i = (— l)"5int(n/2),i for n = 0, 1 ... 7. The quantities A and J 
can be directly determined using the following expressions: 

1 23 

= ^ E Yd (24) 

1 23 y^n 

A. = pE^H^:d-- (25) 

We can then substitute the resulting values of A and J into the linear relation of Eq. 
and test this linear description of our coarse-lattice results for the 24 Cartesian data sets. 
The simplest test of linearity should be Xc '^^ 1^ However, for our 24 sets of seven 
quantities, the resulting Xc/i'^ ' 24) is ~ 15 suggesting this linear description is poor. This 
large value of Xc comes from the linear prediction of the heavy-heavy spin average masses. 
If these are dropped from the calculation of Xc^ we obtain Xc/i^ ' 24) = 1.7, a much more 
acceptable value. Looking more closely, we find the linear prediction for the heavy-heavy 
spin average masses agrees with the calculated value with a fractional discrepancy of 1-2% 
for the 24 data sets. This is certainly a reasonable accuracy given the systematic errors in 
determining these masses from our lattice calculation. However, since the statistical error 
on these quantities, which is used in our definition of Xc^ is of the order of 0.1-0.2% we 
should expect these large Xc values. Thus, we conclude that the linear description of the 
coarse lattice results is satisfactory. 

Using these results for J and A, Eq. [221and the procedure outlined in the previous section 
to determine the error we can go on to find the coarse lattice parameters which describe the 
fine lattice results: 

= {mo, cb, ce, C} = {-0.018(100), 1.648(227), 0.957(904), 1.038(23)}, (26) 

where the superscript T indicates the transpose of the column vector Xq. The results for 
mo and ( are reasonably accurate. Note the relative error in mo should not be determined 
by comparing to the central value for mo which is shifted by the additive renormalization 
implied by merit to be close to zero. Rather, one should recognize that this error in mo 
corresponds to a 4% relative error in m^^. However, the errors on cb and especially ce are 
unacceptably large. 
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In order to better understand these large errors, we now examine the matrix J^-J. This 
matrix is closely related to the matrix J^W^^J which is inverted in Eq. |221to obtain the 
coarse-lattice parameters. While the characteristics of the matrix J^W^^J are entirely 
similar to those of J^-J we found it more natural to focus on the simpler matrix J^-J whose 
definition does not depend on a somewhat ad hoc choice for the correlation matrix WV- 

The eigenvalues of the matrix ■ J are 

{9.55(15), 1.39(10), 0.000138(21), 0.000037(12)} (27) 

with corresponding eigenvectors 

{ 0.832(4), -0.1099(6), -0.1079(7), 0.532(6) }, 



(28) 



{ -0.522(7), 0.062(6), 0.085(3), 0.846(4) }, 
{ 0.181(7), 0.81(7), 0.56(10), -0.003(6) }, 
{ 0.041(23), -0.57(10), 0.82(7), -0.0156(29) } 

Here the eigenvectors reading top to bottom correspond the eigenvalues in Eq. |2ZI reading 
left to right. The eigenvalues span a range of more than five orders of magnitude and 
dramatically decrease between the eigenvectors dominated by the mo or ( directions and 
those aligned with cb or ce- The smallest eigenvalue corresponds to an eigenvector that 
has a large component in the ce direction which leads to large error in the ce coefficient. 
(Recall that the components of the eigenvectors displayed in Eq. EE are arranged in the 
order {mo, cb, Ce, C}- 

Given the range of quantities measured and the precision of the results, we were surprised 
that Ce remains to a large degree undetermined. Of course, this is precisely the result that 
would be obtained if we were working with a redundant set of parameters. Thus, we went 
back and looked carefully at the arguments which determined this set of "independent" 
parameters and discovered an additional field transformation that permits c^; to be trans- 
formed into cr. This result is valid to all orders in ma and up to errors of order (ap)'^. This 

n 

theoretical analysis is presented in the companion paper |15| . 

Here, we will exploit this substantial simplification and use only the three parameters mo, 
( and cp = Cb = Ce- As is shown in the next section, within this restricted parameter space, 
the problem of determining mo, C ^md cp from given values for our seven measured quantities 
is well-posed and accurate results for these three parameters can be easily obtained making 
our proposed step-scaling, matching procedure quite practical. 
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B. Three-Parameter Action 

We will now exploit this simplification from four action parameters to three and determine 
those three parameters which give coarse lattice results agreeing with those found on the fine 
lattice. Specifically, we will use the action in Eq. [T]but fix = = 1 and ce = cb = cp 
and study the dependence of the seven spectral quantities making up the vector Y in Eq. 
on the three parameters mo, cp and C making up the vector: 



(29) 



Cp 

V c J 

As is shown in Ref. a proper, mass-dependent choice for three parameters will yield 
on-shell quantities which are accurate to arbitrary order in (ma)" with errors no larger than 
{apf. 

How does this affect our analysis? We could, of course, disregard all of our four-parameter 
runs and collect an entirely new set of data with the restriction cb = ce- Instead we will 
exploit the approximate linearity of much of our four parameter data and interpolate to 
obtain what we expect to be a good approximation to the results we would obtain had we 
chosen cp = ce- 

Thus, we set cp = cp and exphcitly subtract the deviation that results from ce ^ cp 
using the matrix of derivatives J determined in the four-parameter analysis above. Such 
an expansion in — c^; should be especially safe given the very weak dependence on this 
difference that we have seen. Hence the coarse lattice masses to be used in this three- 
parameter analysis are obtained from: 

'" = y," + J.,2(cS-c^). (30) 

The action parameters corresponding to each of these data sets are Xo^"*'" = , X^"*'" = 
and X2^'*'"' = Xg . The resulting "three-parameter" data sets with 1 < n < 66 can then be 
analyzed in precisely the same fashion as was done for the case of four parameters, following 
the steps taken in Eqs. UHl through |221 

Again we use as the center point that data set giving results closest to the results from 
the fine lattice, which is (X^^).!^)^ ^ {o.0328, 1.511, 1.036} from set #14, the "Cartesian" 
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data sets 43-66, and obtain 

{X^^Y = {mo, cp, C} = {0.037(26), 1.50(9), 1.029(14)}. (31) 

The errors quoted here are statistical and obtained as described in the beginning of this 
section by combining in quadrature the errors coming from the determination of the fine 
lattice masses and the statistical uncertainties in determining the coarse lattice parameters 
which reproduce those fine lattice results. 

Note that mo is relatively small (close to zero) as a reflection of mcriticai for Wilson-type 
fermions lying close to mcharm for our lattice spacing. The significance of the error in mo 
can be estimated from times the error in mo from the average coarse data, giving a 4% 
effect of the error in mo on the resulting heavy-heavy, spin-averaged mass. 

These better defined results for the case of the three-parameter action demonstrate that 
the singularity in the matrix that must be inverted to solve for these heavy quark parameters 
has disappeared. For completeness we list the eigenvalues and eigenvectors of the 3x3 matrix 
(j(3)^Tj(3) contrasted with the singular results found for the four-parameter case in 

Eqs.EZIand|2a 

{9.77(15), 1.41(10), 0.00026(4)}, (32) 
with corresponding eigenvectors 

{ 0.824(4), -0.2157(11), 0.524(6) }, 

{ -0.504(8), 0.142(7), 0.852(4) }, (33) 
{ 0.258(4), 0.9661(11), -0.008(7) }. 

A comparison of Eqs.|2Zland|2Hl"with Eqs. 122 and 1221 shows that the first two large eigenvalues 
and corresponding eigenvectors are changed very little by the reduction from four to three 
parameters. 

Next we would like to examine the contribution to systematic error due to ignoring the 
quadratic terms in our analysis. Using our 24 Cartesian data sets, we can calculate the both 
the first (J- matrix) and second derivatives (a quadratic matrix Q) directly, without using a 
fitting procedure. The resulting simple Taylor expansions around the center point are: 

yn ^ y(3),14 ^ j(3).(^(3),n _ ^(3),14) ^ l(^(3),n _ ^(3),14) . g . (^(3),n _ ^(3),14) ^34) 
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where Q is the 3x3 tensor of second-derivatives and n runs from 43 to 66 (including only the 
Cartesian data sets). We can now estimate how much our resulting parameters X depend 
on the quadratic terms and get a reasonable estimate of the systematic error. 

Using this quadratic approximation, we determine the best-fit, coarse lattice parameters 
Xq""^ by minimizing 

^l^^ = (y^ - r(3)A4 _ J(3).(X(3) _ X(3).14) (35) 
. (Y^ - F(3),14 _ j(3) . (^(3) _ ^(3),14) 

-l/2(xf -X(3).i4)^.Q(3).(xf -X(3).i4)). 

The resuh is (xf V = {mo, cp, (} = {0.034(8), 1.50(3), 1.035(5)}, now including the effects 
of quadratic terms. Comparing these numbers with those in Eq. |^ from the linear ap- 
proximation one sees that the quadratic contributions to the results are buried in statistical 
noise. Therefore, we will not include contributions to the possible systematic errors coming 
the neglect of these quadratic terms in the analysis. 

The systematic errors enter as: (a) We use (ma)^ = 0.2^ or 4% as an estimate of the heavy 
quark discretization errors from domain wall fermion calculation on the fine lattice, (b) The 
remaining RHQ heavy quark discretization effect on the coarse lattice are given by (ap)^ = 
(as(/i = lGeV)maY ~ 0.004. (c) Finally we estimate 1.3% as the systematic error arising 
from the matching of the spatial volumes of fine and coarse lattices. Therefore, adding these 
three systematic errors in quadrature gives our final coefficients: (X^^'*)-^ = {mo,Cp,C} = 
{0.037(26)(13),1.50(9)(6),1.029(14)(40)} where the first error shown is statistical and the 
second systematic. 

In our analysis, we have determined three parameters in the action by requiring that 
seven physical quantities agree between the coarse and fine lattices. Can we match fewer 
physical quantities between the coarse and fine lattice spacing calculations and obtain the 
same result? Table |X] summarizes the results for various choices of the quantities being 
matched. As we can see, all the different choices give consistent values for our three action 
parameters, agreeing within one a. Thus, we have very consistent results for different choices 
of calculated quantities which provides a numerical demonstration of the validity of the heavy 
quark version of the Symanzik improvement program being implemented here. 

Let us focus on two choices of measurements: index "E" using all seven measurements 
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and index "B" using only the results from heavy-heavy data. One might hope that the more 
measurements we include in the analysis, the smaller the resulting errors will be. However, it 
should be recognized that the cost in computer time of making the additional measurements 
involving light quarks is high. As we can see, despite its considerable added cost, the index 
"E" set makes only a small improvement on the statistical errors. It may be more sensible 
to double the number of configurations and focus exclusively on the heavy-heavy system in 
future calculations. 

V. COMPARISONS WITH OTHER APPROACHES 

In this section we compare the parameters tuq, ( and cp determined here for our 1/a = 3.6 
GeV effective heavy quark theory with the similar parameters determined by other methods. 
This serves both as an approximate check of the results determined here and an opportunity 
to compare perturbative and non-perturbative methods. 

A. Determining the bare mass mo 

We first consider the bare mass mo. In most treatments this parameter is related to a 

continuum, "physical" quark mass by a combination of a shift coming from the intrinsic 

chiral symmetry breaking of Wilson fermions and a multiplicative renormalization factor 
7 . 

m{fi) = ZmirriQ - mcnt)/a (36) 

where merit locates the value of mo at which the pion mass vanishes and m(/i) represents a 
continuum quark mass, specified by a renormalization condition imposed at the energy scale 
fi. In the discussion below we will use the MS scheme and /i = 2.0 GeV. We have introduced 
an explicit factor of the inverse lattice spacing in Eq. EHl to give the continuum quark mass 
its proper units. 

We begin by determining the value of m(/i) which corresponds to the mj = 0.2 input 
mass used in our reference, f3 = 6.638 domain wall fermion calculation. For domain wall 
fermions Eq. IHElalso applies but mo should be replaced by m/ and —merit by mres, a measure 
of residual domain wall fermion chiral symmetry breaking that is sufficiently small that it 
will be neglected here. While a quenched, (3 = 6.638, domain wall fermion value for 



26 



is not known, the value Zj^ ~ 1.59 obtained at /3 = 6.0 '21'] may not be too far off. (The 
results presented in Ref j4l| can be used to compare evaluated with the DBW2 action 
at two very different lattice spacings, 1/a = 1.3 GeV and 1/a = 2.0 where a change of less 
than 3% is seen.) Thus, we will assume the calculations described in this paper correspond 
to /i^^(/i = 2.0GeV) = 1.72 GeV. (This large value suggests that our choice for on the 
fine lattice may be somewhat larger than is appropriate for the charm quark mass.) 

To relate this result to the value of itlq expected in our (3 = 6.351 calculation we next 
determine merit. The critical quark mass can be estimated using both perturbative and 
non-perturbative methods. The two-loop, perturbative value for merit for the Wilson gauge 
and clover fermion actions has been obtained in Ref. js^: 

merit = g'^^'^ + g'^^'^ (37) 

/V2 _ 1 

= (-0.1628571 

+ 0.04348303CSW + 0.0180958c^w) (38) 

S(2) = (iV^2 _ 



1 Nf ^ 

-0.017537+ —0.016567+ -f 0.00118618 



+ (0.002601 - ^0.0005597 - ^0.0005459 ) cgw 
V Nc 



1 N 



-0.0001556 + —0.0026226 + -/o.0013652 cl^ 

iV 2 Nc J 



1 Nf \ 

-0.00016315 + ^0.00015803 - ^0.00069225 ] c: 



,3 

sw 



1 Nf 
" /, 



+ -0.000017219 + —0.000042829 - -f 0.000198100 c. 



iV2 N, 



c 



c / 



(39) 



where the number of fermion flavors Nf = in the quenched approximation and the number 
of colors Nc = 3. We can obtain the coefficient, cgw, of the clover term from the non- 
perturbative result of Ref. : 

1 - 0.656^2 _ o.i52c,4 _ 0.054C/6 
= 1 - 0.922,2 ' (40) 

where g is the lattice coupling constant and P = 6/g'^. This gives csw = 1.544 on a /3 = 6.351 
lattice and merit = —0.219. 

An alternative way of computing merit is to take the non-perturbative, the ALPHA col- 
laboration measurement of Kerit (for example from Table 1 in Ref. Q]) and parameterize it 
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function of coupling constant: 



0.130287 - 0.239546C/2 + 0.111829c/^ 
~ 1 - 1.84915C/2 + 0.868181^4 ^ ) 



for 6.0 < P < 7.4 and use merit = 4. This gives merit = —0.317. We will adopt this 

^^crit 

latter, non-perturbative value as being more accurate. 

Finally, in order to invert Eq. IHSl to obtain the expected value of mo which can be 
compared with our results we require the appropriate factor Zm for our rather fine f3 = 
6.351 lattice and relatively large, m^^{fi = 2.0GeV) = 1.72 GeV quark mass. For this 
comparison, we can avoid the extra translation to and from the MS scheme by directly 
comparing quantities calculated in the RI scheme at = 2.0 GeV. From Tables I and 
II of Ref. j2l| we determine Z^^(DWF) = 1.81. We will use a similar non-perturbative 



value ^^^(SW) = = 1.82 extracted from Table 1 of Ref. |43|. This value is only 

approximate for our situation since it was obtained for light quark masses and on a coarser, 
P = 6.2, lattice. Using these values we obtain mo = .^^^(DWF)/Z^'^(SW) ■ m/ + merit = 
—0.018 in units of 1/a = 3.6 GeV. Since the light quark value for Z^^{SW) chosen in this 
estimate will have 0{ma) errors, we expect systematic errors of size ~ 0.1, and should view 
the agreement between this mo = —0.018 estimate and the mo = 0.036(34) result in Eq. |^ 
as very reasonable. 

B. One-loop RHQ coefficients 

We now compare our non-perturbative result for the remainin g pa rameters ( and cp with 
the one-loop perturbative calculations carried out by M. Nobes |26| for the closely related 
quantities, (, cb and ce appearing in the Fermilab action. 

These one-loop coefficients of Fermilab action were calculated using automated per- 
turbation theory techniques from the scattering of a quark off of a background chromo- 
magnetic (electric) field j2a]. The calculations are done on the lattice and in the continuum 
and the comparison used to determine the lattice parameters. 

The analytic tree-level coefficients (after being translated into our notation for the action) 
are: 



_ ' mo(2 + mo) y mo(2 + mp) mo(2 + mo) 
^ ~\ 1^4(1 + mo); 21n(l + mo) 4(1 + mo) 
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e^l = CM ( ^^"^''^ + + "^°(^ + "^°) ^ (42) 

Next, the one-loop result for is given by the formula: 

- ri , .^[1] A (C'°')^ + C'°'sinh(ln(l + mo)) 
C - l^ + ^o^M.J C[o] + sinh(ln(l + mo)) ' ^ ^ 

We use this formula and the numerical one-loop results from Tables 6.1 and 6.2 in Ref. [3| 
and perform an error- weighted fit to the three functions of interest, ^(^0*^), CB(moa) and 
c^(moa) with expressions of the form 

1 + ^ ' 



where X represents cs, ce and C, while the a's and (i's are listed in Table IXII This 
fit implies that at tjiq = 0.036, the coefficients are cb = 1.261, ce = 1.246 and 
( = 1.003, ~ lAa lower than our non-perturbatively determined coefficients: (Xq^^)'^ = 
{0.037(26)(13), 1.50(9)(6), 1.029(14)(40)}. (Since the results of Nobes have cb ~ ce we can 
directly compare the coefficients in his 4-parameter and our 3-parameter lattice action.) 

To see directly the effects of the differences between these perturbative and non- 
perturbative coefficients, we should compare the resulting spectra. Although we did not 
use these one-loop numbers in a spectrum calculation, we can use our linear description of 
the dependence of the spectra on the action parameters (the coefficients J and A of Eq. ITH)) 
to get a good idea of what the resulting masses would be were we to use these one-loop 
coefficients. We summarize the results in Table IXITl These are reasonably close to our non- 
perturbative results with the largest discrepancy being the two hyperfine splittings which 
are 25% smaller when determined from the one-loop coefficients. 

There is a second, extensive perturbative calculation of the one-loop, tadpole improved 
RHQ lattice action by the Tsukuba group [3]. However, because the Tsukuba action uses 
five parameters with 7^ C, we cannot make a direct comparison. While continuum field 
transformations can be employed on the continuum effective Lagrangian to prove that these 
5-parameter and 3-parameters descriptions should lead to the same continuum physics up to 
discrepancies of order (pa)^, these transformations are not defined for the lattice variables 
and cannot be used to relate the one-loop coefficients of the Tsukuba action given in Ref. [3| 
and those determined here. 
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VI. SUMMARY AND OUTLOOK 

In this work, we have demonstrated that it is practical to determine the coefficients of 
the relativistic heavy quark action non-perturbatively through a finite-volume, step-scaling 
technique. This has been done by matching various heavy-heavy and heavy-light mass 
spectrum calculations on two quenched lattices. The domain wall fermion action is used on 
fine lattice, where ma is relatively small, while for the coarse lattice an improved relativistic 
heavy quark action is used. By comparing the finite- volume predictions of these two theories 
we can then determine the coefficients of the heavy quark action. In order to simplify 
the analysis, we assumed a linear relation between the parameters appearing in the heavy 
quark action and the resulting mass spectrum. The coefficients in this linear relation were 
determined by computing the coarse-lattice mass spectrum for a number of choices for the 
RHQ action. We could then use this linear relation to precisely determine those heavy quark 
parameters which would yield the masses implied by the fine lattice calculations. 

We initially applied this matching technique to the four-parameter version of the heavy 
quark action originally proposed in Ref. However, for this case, the system of lin- 

ear equations that must be solved was singular within statistical errors and the resulting 
parameters, especially the coefficients cb and ce very poorly determined. This lead us 
to search for possible redundancy in the four-parameter action and recognize, as is dis- 
cussed in detail in a companion paper that a further field transformation was avail- 
able that could be used to set ce = cb, reducing the number of independent param- 
eters to three. With this restriction the problem of determining the relativistic heavy 
quark action is well-posed and the coefficients can be accurately determined. Our result 
for the bare mass, clover term and asymmetry between the space and time derivatives is 
{mo,cp,C} = {0.037(26)(13), 1.50(9)(6), 1.029(14)(40)}, where the first error is statistical 
and the second systematic, excluding those coming from the quenched approximation. Fi- 
nally, we included a quadratic term in the dependence of our measured masses on the action 
parameters and obtained a result consistent with the linear expansion. 

We can easily decrease the statistical error by increasing the number of configurations 
(here 100 were used) and reduce the systematic error by starting with a finer lattice for 
the domain wall fermion calculation. Our use of the quenched calculation is intended to 
provide a computationally inexpensive study of the matching procedure. The next step is 
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a determination of the coefficients in this relativistic heavy quark action, appropriate for 
charm physics in full QCD. As discussed in Sec. lUl we can perform the same finite- volume, 
step-scaling procedure on 2+1 flavor dynamical lattices. Since the long- and short- distance 
physics can be treated separately, we can substantially reduce the computational cost of 
such full QCD step scaling by using heavier light quark sea masses in the earlier stage of 
matching, as long as msf^a/ ^qcd are equal for each pair of systems being matched. Such a 
calculation should be practical on presently available computers. 
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TABLE I: Comparison of the conventions and/or values used to specify the six terms in the 
improved lattice action of Eq. ^ The top row identifies terms that appear in Eq. ^ The next row 
lists our choice for the coefficient of each term and the next two rows specify that same coefficient 
written in the notation of the Fermilab and Tsukuba 13] papers. 



Action 




7B 




-0? 


irr- F- ■ 




This paper 


1 


c 


1 


rs 


cb 


ce 


Fermilab 


1 


c 


1 


rsC 


cbC 


ceC 


Tsukuba 


1 


V 


1 


rs 


CB 


CE 



TABLE IL The specific choice of parameters for the two sets of lattice configurations analyzed in 
this paper. 

I II 



Volume X 48 16^ X 32 
1/a 5.4 GeV 3.6 GeV 
L 0.9 fm 0.9 fm 

/3 6.638 6.351 

ma 0.2 0.3 

Action DWF RHQ 
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TABLE III: The choice of parameters for the two sets of lattice configurations needed for the next 
step in this step-scaling program. Calculations with these parameters are now underway but are 
not described in this paper. 

Ill IV 

Volume X 48 16^ x 32 
1/a 3.6 GeV 2.4 GeV 
L 1.33 fm 1.33 fm 

(3 6.351 6.074 

ma 0.3 0.45 

Action RHQ RHQ 



TABLE IV: Meson states created by local operators of the form il^Vil}, labelled in spectroscopic 
notation. 

r ''s+^Lj 

75 0-+ 

li ^s, 1— 

1 ^Po 0++ 

757i ^Pi 1++ 



TABLE V: Common parameters for each of the coarse and fine data sets. For both data sets 
L = 0.9 fm while for the domain wall fermion action we use = 12 and M5 = 1.5. Here "TBD" 
indicates a value to be determined in the matching proceedure being developed here. 

Label (3 \ Sl arriL Sh amn a~^{SQ, pot.) 



Fine 6.638 24^ x 48 DWF 0.02 DWF 0.2 5.4 GeV 
Coarse 6.351 16^ x 32 DWF 0.03 RHQ TBD 3.6 GeV 
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TABLE VI: Mass spectrum measured on the fine lattice in units of a 

mps my nil/ 1712 rn ay ms 

light-light 0.175(3) 0.233(5) _ _ _ 

heavy-light 0.467(2) 0.485(3) _ _ _ 

heavy-heavy 0.716(1) 0.728(1) 1.02(2) 0.810(5) 0.799(4) 



TABLE VII: Light-light hadron spectrum measured on the coarse lattice in units of a and units 
of 3/2a to compare with the fine lattice results in Table IVll 

units JTips my 



light-light 1/a 0.259(6) 0.328(10) 
light-light 3/2a 0.173(4) 0.219(7) 
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TABLE VIII: Parameters used on the coarse lattice. 



Set # 
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0. 


1 ^^20f> 

J- .tj'j ^yjyj 


1.45769 


1.01281 
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.01747 
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1 00063 


35 
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1.00674 
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37 
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TABLE IX: Mass spectrum measured on the coarse lattice in units of (where "sa" is spin 
averaged, "hs" is hyperfine spHtting, "soa" is spin-orbit averaged, "sos" is spin-orbit sphtting; "hi" 
is heavy-light and "hh" is heavy-heavy). 



Set # 



1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 



1.0137(17; 
1.1314(15 
1.0872(16; 
1.1755(16 
0.9776(16 
1.0724(17; 
1.0840(15; 
1.0622(15 
1.1177(15 
1.0949(14 
1.0255(15; 
1.1148(15; 
1.0981(15 
1.0937(15 
0.9384(19 
0.9637(18 
0.9822(18 
0.9935(18 
1.0115(18 
1.1019(16 
1.1187(16 
1.1353(16 
1.0809(16 
1.0521(16 
0.9985(17 
0.9143(19 
0.7896(23 
0.606(3) 
0.8593(22 
0.9229(20 
0.9461(19 
0.9573(19 
0.9693(18 



0.0174(6 
0.0152(5 
0.0160(5 
0.0153(5 
0.0169(5 
0.0171(5 
0.0170(5 
0.0176(5 
0.0165(5 
0.0161(5 
0.0174(5 
0.0157(5 
0.0160(5 
0.0162(5 
0.0204(7 
0.0199(6 
0.0195(6 
0.0193(6 
0.0190(6 
0.0171(5 
0.0169(5 
0.0167(5 
0.0164(5 
0.0175(6 
0.0195(6 
0.0229(8 
0.0295(11 
0.0469(28 
0.0235(8 
0.0212(7; 
0.0204(7; 
0.0201(7 
0.0197(6 



0.679(5 
0.742(4; 

0.718(5 
0.762(5 
0.664(5 
0.709(5 
0.716(4; 
0.704(4; 

0.734(4; 
0.723(4; 

0.686(4 
0.733(4; 

0.724(4; 
0.721(4; 
0.638(5 
0.652(5 
0.661(5 
0.667(5 
0.677(5 
0.725(5 
0.734(5 
0.743(5 
0.714(5 
0.699(5 
0.670(5 
0.624(5 
0.557(5 
0.460(6 
0.596(5 
0.630(5 
0.642(5 
0.648(5 
0.655(5 



0.0254(17 
0.0221(15 

0.0233(15 
0.0228(15 
0.0241(16 
0.0248(16 
0.0246(16 
0.0254(16 
0.0238(15 
0.0234(15 
0.0254(16 
0.0232(15 
0.0236(15 
0.0239(15 
0.0287(19 
0.0280(18 
0.0275(18 
0.0273(18 
0.0268(17 
0.0242(16 
0.0239(16 
0.0235(15 
0.0240(16 
0.0255(16 
0.0279(18 
0.0318(20 
0.0384(25 
0.052(4) 
0.0316(21 
0.0293(19 
0.0286(19 
0.0282(18 
0.0279(18 



0.020(5; 
0.017(4 

0.018(4; 
0.017(4; 
0.020(5 
0.019(5; 
0.018(4 
0.025(4 
0.017(4 
0.017(4 
0.019(4 
0.017(4; 
0.017(4 
0.018(4 
0.024(6 
0.023(5 
0.022(5 
0.022(5 
0.021(5 
0.019(5; 
0.018(4 
0.018(4 
0.019(5 
0.020(5 
0.023(6; 

0.027(7; 

0.035(9 
0.053(18) 
0.026(7 
0.024(6; 
0.023(6 
0.023(6 
0.023(5 



1.135(14; 
1.251(13 
1.207(14 
1.298(13 
1.094(14 
1.193(14 
1.208(12 
1.190(12 
1.241(12 
1.217(12 
1.149(12 
1.238(12 
1.221(12 
1.216(12 
1.060(15 
1.085(15; 
1.103(15 
1.114(14 
1.132(14 
1.221(13; 
1.237(13 
1.254(13 
1.196(14 
1.164(14 
1.107(14; 
1.021(15 
0.895(18 
0.708(24 
0.993(17 
1.050(15 
1.070(15; 
1.079(15 
1.089(15 



0.988(21) 
0.949(18) 
0.966(19) 
1.075(21) 
0.842(16) 
1.021(21) 
1.009(17) 
1.016(18) 
1.002(16) 
0.956(15) 
0.980(17) 
0.992(16) 
0.995(16) 
0.998(16) 
1.016(23) 
1.013(23) 
1.012(22) 
1.011(22) 
1.009(22) 
0.986(20) 
0.985(19) 
0.984(19) 
0.973(19) 
0.988(20) 
1.013(22) 
1.054(27) 
1.12(4) 
1.21(10) 
1.049(28) 
1.027(25) 
1.019(24) 
1.016(23) 
1.011(23) 
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Set # 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 



0.9760(18 
1.2455(15 
1.0711(17; 
1.0736(17; 
1.0696(17; 
1.0700(17 
0.9864(18 
1.1359(16; 
1.0824(17; 
1.1866(15; 
0.8149(23 
1.2271(15 
0.8705(22 
1.3000(14; 
0.9691(19; 
1.3356(14 
1.0164(18; 
1.1473(15 
0.7632(24 
1.1894(15 
0.8218(22 
1.2662(14 
0.9272(19; 
1.3030(14 
0.9764(18 
1.2495(15 

0. 9026(20 
1.0626(17; 
1.1071(17; 
1.0140(18 
1.1453(16 
1.1051(17 

1. ()()5()(l 



0.0195(6 
0.0136(4 
0.0171(5 
0.0170(5 
0.0172(5 
0.0176(5 
0.0198(6 
0.0162(5 
0.0173(5 
0.0157(5; 
0.0245(9 
0.0138(4; 
0.0212(8; 
0.0136(4; 
0.0189(6; 
0.0119(3 
0.0165(5 
0.0160(5 
0.0257(9; 
0.0141(4; 
0.0221(8 
0.0138(4; 
0.0194(6; 

0.0120(3; 

0.0169(6 
0.0137(4; 
0.0200(7 
0.0172(5 
0.0151(5; 
0.0179(6 
0.0148(5 
0.0160(5 



0.658(5 
0.802(4; 
0.708(5 
0.709(5; 
0.707(5 
0.708(5 
0.663(5 
0.743(5 
0.715(5 
0.769(4; 
0.571(5 
0.790(4; 
0.601(5 
0.830(4 
0.654(5 
0.849(4; 
0.679(5 
0.749(4 
0.545(5 
0.772(4; 
0.576(5 
0.813(4; 
0.633(5 
0.833(4 
0.659(5 
0.803(4; 
0.619(5 
0.704(4; 
0.728(4; 
0.678(4; 
0.748(4; 
0.726(4; 
0.70()(1 



0.0277(18 
0.0195(13 
0.0248(16 
0.0248(16 
0.0248(16 
0.0249(16 
0.0279(18 
0.0228(15 
0.0245(16 
0.0227(15 
0.0340(23 
0.0209(14 
0.0309(21 
0.0195(13 
0.0274(18 
0.0179(12 
0.0250(17 
0.0230(15 
0.0354(24 
0.0210(14 
0.0319(21 
0.0196(13 
0.0279(18 
0.0179(12 
0.0254(17 
0.0195(13 
0.0284(19 
0.0242(16 
0.0221(15 
0.0256(17 
0.0212(14 
0.02.30(15 
0.02:«(1."') 



0.023(5 
0.013(4 
0.019(5 
0.019(5 
0.019(5 
0.019(5 
0.022(5 
0.017(4 
0.019(5 
0.017(4; 
0.030(7 
0.015(4 
0.027(6 
0.014(4 
0.022(5 
0.012(3 
0.020(5 
0.018(4 
0.032(8 
0.016(4 
0.028(7 
0.014(4 
0.023(6 
0.013(3 
0.020(5 
0.014(4 
0.023(6 
0.018(5 
0.016(4; 
0.020(5 
0.015(4 
0.017(4 
0.017(1 



1.095(14 
1.372(13 
1.192(14 
1.194(14 
1.191(14 
1.194(14 
1.108(15 
1.262(13 
1.205(14 
1.303(13 
0.937(17 
1.344(13 
0.993(16 
1.420(12 
1.095(15 
1.456(12 
1.143(14 
1.263(13 
0.884(17; 
1.305(13 
0.943(16 
1.385(13 
1.052(15 
1.422(12 
1.102(15 
1.367(12 
1.027(15 
1.184(14 
1.228(13 
1.134(14 
1.267(13 
1.226(13 
1.18()(l:-! 



m\/m2 
1.009(22) 
0.976(18) 
1.020(21) 
1.022(21) 
1.019(21) 
1.013(21) 
1.036(23) 
0.997(20) 
1.007(21) 
1.019(20) 

1.10(3) 
1.012(19) 
1.090(28) 
0.986(18) 
1.052(24) 
0.979(17) 
1.045(22) 
0.973(19) 

1.06(3) 
0.967(18) 
1.045(28) 
0.940(17) 
1.004(23) 
0.934(17) 
0.997(21) 
0.978(18) 
1.052(24) 
1.016(20) 
1.009(19) 
1.034(21) 
0.995(19) 
1.0.36(20) 
0.<)89(19) 
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TABLE X: The resulting coarse-lattice paramters obtained by matching various combinations of 
physical quantities between the coarse and fine lattices. Different choices of for the quanties to be 
matched give action parameters consistent with each other within one cr, showing the consistency 
of this heavy quark improvement program. 

Measurement Index Data used X^^ 

A m^^, m'^^, m^^^, mi/m2 { 0.07(4), 1.67(13), 1.030(14)} 

B "A" +m^^^ { 0.04(3), 1.56(10), 1.034(12)} 

C m^i*, m'^^, m^i, rn^^l, mi/m2 { 0.06(4), 1.62(13), 1.032(14)} 

D "D" +m^^, { 0.04(3), 1.56(10), 1.034(12)} 

E ah { 0.03(3), 1.53(10), 1.035(12)} 



TABLE XI: Parametrization of the one-loop coefficients of the Fermilab action using Eq. 1441 



ao ai 


02 


as 


di d2 dz 


C'^] 0.00029923 0.00124977 


0.163759 


0.0258287 


5.10243 1.65713 0.00633212 


c|' 0.270419 0.431474 


0.162718 


0.00212438 


1.87436 0.319194 0.00619183 


4' 0.271519 0.0122322 


-0.000039117 





0.0565955 
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TABLE XII: The expected coarse lattice results for various choices of coefficients, Xq, in the 
heavy quark effective action. Here the linear approximation of Eq. 1181 with coefficients A and J 
determined from data sets #43-66, is being used to predict the corresponding physical masses. 



Parameter 










^hh 

'""SOS 


^hh 
'"soa 


mi/m2 



X^^-"' 1.0854(16) 0.0165(5) 0.716(4) 0.0239(16) 0.019(4) 1.206(13) 1.000(20) 
XFermiiab 1.1700(15) 0.0124(4) 0.763(4) 0.0194(13) 0.014(4) 1.292(13) 0.943(16) 
scaled Yjr 1.0881(16) 0.0175(6) 0.725(3) 0.0267(13) 0.017(3) 1.211(8) 1.002(23) 
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FIG. 1: The sequence of lattice sizes and lattice spacings used to determine the coarse-lattice, 
heavy quark parameters through a step-scaling technique beginning with a comparison with an 
0(a)-improved light quark calculation. The matching between the top two lattice spacings is the 
calculation described in this paper. 
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FIG. 2: The auto-correlation function for five different heavy-heavy meson propagators evaluated 
at a source-sink separation of twelve lattice units with /? = 6.638 and a 24^ x 48 space-time 
volume. In the top graph the propagators were calculated on every 5th configuration while in the 
bottom graph the propagator measurements were separated by 50 sweeps. This suggests that our 
separation of 10,000 sweeps between measurements ensures that they will be uncorrelated. 
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FIG. 3: The static quark potential calculated on both the coarse and fine lattices. The static 
quark potential values computed on the /5 = 6.351 lattice are shown as circles. The squares mark 
the static quark potential from the /3 = 6.638 lattice scaled by the fitted ratio of lattice spacings 
A = 1.51 and shifted by a constant. The agreement between these two different sets of points gives 
good evidence that the ratio of lattice spacings between these two /3 values is the desired 3/2. 
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FIG. 4: The dependence on the fifth dimension, s, of the space-time sum of the modulus squared 
of the nineteen lowest-lying eigenvectors of the hermitian domain-wall Dirac operator: Ylix l*a;,sP- 
This is shown for a Wilson quenched gauge configuration with = 6.638, V = 16^, Lg = 12, 
Ms = 1.8, and m/ G {0.22,0.27,0.37,0.47}. The unphysical, propagating states are seen only for 
mf> 0.37. 
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FIG. 5: A comparison of effective mass plateaus between point and smeared sources in the heavy- 
heavy sector. Triangles denote the point-point source, and pentagons denote the point-smeared 
heavy meson correlators. Clearly, the plateaus have been improved by the smearing. However, a 
double-cosh fit to the two distinct wavefunction sources might help us determine the ground state 
energy even more accurately. 
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FIG. 6: The distribution in the 3-parameter space {rriQ, cb, C) of the 24-set "Cartesian" data. 
The center circular point is set ^^14, and the points around it are sets#43-66. The starred point 
represents the final matching coefficients determined in Sec. II V Bl 



